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Abstract 

We consider supergravity with a gravitino lightest supersymmetric particle. The next-to-lightest 
supersymmetric particle (NLSP) decays to the gravitino with lifetime naturally in the range 10"^ — 
10^ s. However, cosmological constraints exclude lifetimes at the upper end of this range and 
disfavor neutralinos as NLSPs, leaving charged sleptons with lifetimes below a year as the natural 
NLSP candidates. Decays to gravitinos may therefore be observed by trapping slepton NLSPs in 
water tanks placed outside Large Hadron Collider (LHC) and International Linear Collider (ILC) 
detectors and draining these tanks periodically to underground reservoirs where slepton decays 
may be observed in quiet environments. We consider 0.1, 1, and 10 kton traps and optimize their 
shape and placement. We find that the LHC may trap tens to thousands of sleptons per year. 
At the ILC, these results may be improved by an order of magnitude in some cases by tuning the 
beam energy to produce slow sleptons. Precision studies of slepton decays are therefore possible 
and will provide direct observations of gravitational eff'ects at colliders; percent level measurements 
of the gravitino mass and Newton's constant; precise determinations of the gravitino's contribution 
to dark matter and supersymmetry breaking's contribution to dark energy; quantitative tests of 
supergravity relations; and laboratory studies of Big Bang nucleosynthesis and cosmic microwave 
background phenomena. 

PACS numbers: 04.65. +e, 12.60.Jv, 26.35. +c, 98.80.Es 
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I. INTRODUCTION 



Weak-scale supersymmetry remains a beautiful framework for resolving the problems of 
electroweak symmetry breaking and dark matter and its discovery is among the most 
eagerly anticipated events in particle physics. Opportunities for supersymmetry discoveries 
and studies at colliders depend largely on which superpartner is the lightest supersymmetric 
particle (LSP). In simple supergravity models, supersymmetry is transmitted to standard 
model superpartners through gravitational interactions, and supersymmetry is broken at a 
high scale. The mass of the gravitino G is 

F 

and the masses of standard model superpartners are 

where F ~ (10^^ GeV)^ is the supersymmetry breaking scale squared, and = 
(SttCtv)"^/^ ^ 2.4 X 10^^ GeV is the reduced Planck scale. The precise ordering of masses 
depends on unknown, presumably C(l), constants in Eq. Most supergravity studies 
assume that the LSP is a standard model superpartner, such as a slepton or neutralino. 

Recently attention has turned to the other logical possibility, namely, that the gravitino 
is the LSP jli 0, 0, 0, 0, 0, S, 0, • In supergravity where supersymmetry breaking 

is mediated by gravity, the gravitino has a mass ~ 100 GeV if the superpartner mass 
scale is m ~ 100 GeV. The gravitino couplings are also suppressed by M^,. The gravitino's 
extremely weak interactions imply that it is irrelevant for most supersymmetric processes. 
However, if the gravitino is the LSP, the next-to-lightest supersymmetric particle (NLSP) 
decays to its standard model partner and a gravitino. The NLSP is a weak-scale particle 
decaying gravitationally and so has a natural lifetime of 

^ ~ 10^ - 10^ s , (3) 

as will be discussed more fully in Sec. [H] This lifetime emerges naturally in this simple 
supersymmetric scenario. At the same time, it is outlandishly long by particle physics 
standards It requires a revamping of many aspects of supersymmetric phenomenology 
and cosmology and opens up novel opportunities, including the one discussed here. 

The gravitino LSP scenario is constrained by cosmological and astrophysical data. The 
gravitino is a stable superweakly massive particle (superWIMP) and forms dark matter. Its 
production during reheating and by NLSP decays is therefore constrained by measurements 
of the non-baryonic cold dark matter density. NLSP decays also deposit electromagnetic 0, 
Esl . [l6l [Tt^ and hadronic [3, energy into the universe well after Big Bang nucleosynthesis 



(BBN), and so may ruin the successful predictions of standard BBN. These decays may also 
distort the cosmic microwave background (CMB) from its observed Planckian spectrum. 
Last, photons produced in NLSP decays are subject to bounds on the diffuse photon flux. 

The impact of these constraints on the gravitino LSP scenario have been considered 
in detail. In addition to the leading two-body NLSP decays to the gravitino 0, IE B 
three-body NLSP decays must also be considered when they are the leading contribution to 
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hadronic cascades The result is that the gravitino LSP scenario is not excluded and, 

in fact, all constraints may be satisfied for natural weak-scale NLSP and gravitino masses. 

Not all possibilities are allowed, however, and two results are particularly worth noting. 
First, neutralino NLSPs are highly disfavored Neutralinos typically have two-body 

decays x ~^ ZG qqG. The resulting hadronic cascades destroy BBN successes, and 
exclude this scenario unless such decays are highly suppressed. Kinematic suppression is 
not viable, however — if — ttlq < mz, the decay x ~^ iG takes place so late that it 
violates bounds on electromagnetic cascades. Neutralino NLSPs are therefore allowed only 
when the two-body decays to Z bosons are suppressed dynamically, as when the neutralino 
is photino-like, a possibility that is not well-motivated by high energy frameworks. Slepton 
and sneutrino NLSPs also produce hadronic energy when they decay, but this occurs only 
through three-body decays. These have been analyzed and found to be safe [2,, As a 
result, the most natural NLSP candidates are sleptons, particularly the right-handed stau. 

Second, cosmological constraints exclude the upper range of lifetimes in Eq. (jS)) 0, 01 ■ 
Very late decays occur in a cold universe where decay products are not effectively thermalized 
and so are especially dangerous. For typical thermal relic NLSP abundances, the CMB and 
BBN constraints therefore provide an upper bound on NLSP lifetimes, roughly excluding 
those above a year. 

In passing, we note that the scenario outlined above has a number of other motivations. 
One such motivation is from BBN. Late NLSP decays not only pass all BBN constraints, 
they may even resolve the leading BBN anomaly by destroying ^Li to bring the predicted 
abundance in line with the low values favored by observations 0, I3- To resolve the ^Li 
anomaly, the preferred NLSP lifetime is ~ 3 x 10^ s j^, that is, about a month. A second 
motivation follows from considerations of leptogenesis ji^. Gravitinos may be produced 
during reheating. If the gravitino is not the LSP, its late decays are dangerous to BBN, 
and require reheating temperatures Trh ~ 10^ GeV to 10^ GeV [l9|, in conflict with the 



requirement Trh ~ 3 x 10^ GeV of thermal leptogenesis 21[. In contrast, in the gravitino 
LSP scenario, the gravitino does not decay, and the reheat temperature is bounded only by 
the overclosure constraint on the gravitino density. For ttiq ~ 100 GeV, reheat temperatures 
as high as ~ 10^° GeV are allowed |23|], consistent with thermal leptogenesis |24 i25t . 
Additional connections between leptogenesis and gravitino LSPs are discussed in Ref. |26l |. 

Given all of these motivations, we investigate here the collider implications of a gravitino 
LSP with a charged slepton NLSP with lifetime under (but not much under) a year (27l |. 
In particular, we investigate the possibility of trapping sleptons in material placed just 
outside Large Hadron Collider (LHC) or International Linear Collider (ILC) detectors. This 
material may then be moved to some quiet location so that slepton decays may be observed 
in a relatively background- free environment. Although these objectives may be realized in 
many ways, we study here the particularly simple possibility of trapping sleptons in water 
tanks which may be drained periodically to underground reservoirs where the slepton decays 
may be observed. 

In Sec. ini we discuss the relevant properties of sleptons in the gravitino LSP scenario. 
In Sec. mil we discuss our procedure for maximizing the number of sleptons trapped given 
a flxed volume of water. This is applied to the cases of the LHC and ILC in Sees. II VI and 
IVl respectively. At the LHC, we flnd that tens to thousands of sleptons may be trapped 
each year, depending on the overall mass scale of supersymmetry. At the ILC, typically 
far fewer sleptons are produced. However, by controlling the beam energy, the velocity of 
produced sleptons may be tuned to some low value, allowing a large fraction of sleptons to 
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be trapped. By exploiting this feature, we find tliat an order of magnitude more sleptons 
may be trapped at the ILC than at the LHC These results imply that percent level studies 
of slepton decays may be possible. Such studies will have fundamental implications for 
supergravity, supersymmetry breaking, dark matter, and dark energy. These implications 
and our conclusions are discussed in Sec. IVll 



II. SLEPTON PROPERTIES IN THE GRAVITINO LSP SCENARIO 
A. Slepton Mass 

In supergravity with a gravitino LSP, the slepton NLSP is expected to have a weak-scale 
mass. Current collider bounds require mj > 99 GeV from null searches for long-lived charged 
tracks at LEP II 

Cosmology brings additional considerations, however. Gravitinos produced in the late de- 
cays of sleptons are superWIMP dark matter. Barring the possibility of entropy production 
after slepton freeze out, the gravitino relic density must therefore satisfy 

n^h^ = "^nfh^ < fioM/i' , (4) 

mi 

where VL^ is the slepton's thermal relic density, and the non-baryonic cold dark matter 
density is constrained to the range 0.094 < ^dm^^ < 0.124. Assuming itlq and mi are 
not too disparate, this provides an upper bound on the slepton and gravitino masses, since 
oc m?. Without special effects, ~ ^dm for m^ ~ 700 - 1000 GeV. If ^ > 0.1, the 
overclosure constraint requires superpartner masses below about 3 TeV. 

On the other hand, a particularly attractive possibility is that gravitino superWIMPs 
are most or even all of the non-baryonic dark matter. Although not a strict requirement, 
one might therefore prefer i^^M- The decay to the gravitino only reduces the relic 

density. Without special effects, then, overclosure requires mi ^ 700 — 1000 GeV. Such 
heavy sleptons will be difficult to explore at the LHC and are kinematically inaccessible in 
the first stage of the ILC. 

Just as in conventional neutralino dark matter scenarios, however, there are supplemen- 
tary mechanisms for gravitino production. One such mechanism is co-annihilation. If a 
neutralino x is just slightly heavier than the slepton, it will freeze out with the slepton and 
later decay to the slepton, adding its thermal relic density to the slepton's. This allows 
supersymmetric models with much lower slepton masses to produce the correct gravitino 
superWIMP dark matter density. For example, in minimal supergravity with Aq = 0, 
tan/3 = 10, and /i > 0, the desired relic density may be achieved near the f LSP-x LSP 
border at M1/2 = 300 GeV, where ^ m^ ^ 120 GeV [29]. Alternatively, gravitinos may 
be produced during reheating. For reheating temperatures Trh ~ 10^ GeV, as might be 
preferred for leptogenesis as discussed above, gravitinos may again be all of the non-baryonic 
dark matter for slepton masses as low as 120 GeV 2|. It is clear that such effects may be 



very important. We will consider a variety of slepton masses below, including those within 
reach of a first stage ILC. We reiterate that all of these considerations depend on the as- 
sumption that gravitinos make up all of the dark matter. From a purely particle physics 
viewpoint, this is clearly optional — some other particle, such as the axion, may be the dark 
matter. In this case, gravitino LSPs and light NLSP sleptons are perfectly possible without 
any additional restrictions. 
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B. Slepton Lifetime 



The width for the decay of a slepton to a gravitino is 
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(5) 



assuming the lepton mass is neghgible. This decay width depends on only the slepton mass, 
the gravitino mass, and the Planck mass. In many supersymmetric decays, dynamics brings 
a dependence on many supersymmetry parameters. In contrast, as decays to the gravitino 
are gravitational, dynamics is determined by masses, and so no additional parameters enter. 
In particular, there is no dependence on left-right mixing or flavor mixing in the slepton 
sector. 



C. Slepton Range in Matter 



Last, it will be crucial to this study to know the range of sleptons in matter. Charged 
particles passing through matter lose energy by emitting radiation and by ionizing atoms. 
At lower energies, ionization dominates the energy loss, while at high energies, radiation is 
the dominant effect. As we will see below, it is unreasonable to expect to stop sleptons with 
momenta much larger than their rest mass. For the present case, then, ionization losses are 
dominant, and radiation is negligible. 

The average energy loss due to ionization is given by the Bethe-Bloch equation. The low 
energy approximation to the Bethe-Bloch equation may be derived [s^ by first calculating 
the classical cross section for a collision with fixed impact parameter and energy loss. One 
then integrates the impact parameter from the Compton wavelength of the free electron as 
seen by the charged particle to a maximal impact parameter where the particle cannot "see" 
the electron in the time that it passes by the atom. A more refined treatment yields the 
Bethe-Bloch equation [sH 
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(6) 



where dE/dx is the energy loss per g cm"^, K = 0.307075 MeV g~^ cm^, and mg is the 
electron mass. The material is characterized by its atomic charge Z in units of e; its average 
nucleon number A; and its mean ionization energy J, which is given for various elements in 
Ref. The incoming particle has mass M, charge z in units of e, velocity /?, and dilation 
factor 7 = {1 - 

The parameter 6 accounts for the fact that incoming particles polarize the surrounding 
medium. At high energies, this effect may be included by setting 



i = e(B - Bo) 



(7) 



where ujp is the plasma frequency, and £"0 is the energy at which the effect of polarization 
is significant. This correction is typically significant for jS'j ^ 10. We include this effect, 
although, as we will see, sleptons that may be trapped in a reasonably sized detector have 
P'j < 1, and so this effect is also irrelevant for the final results of this study. 
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TABLE I: Range parameters for lead and water, assuming a 219 GeV slepton. 



Material 


6M (MeV) 


/(eV) 


Eo (TeV) 


hiOp (eV) 


Lead 


110 


820 


4.4 


61 


Water 


220 


75 


1.3 


21 



At low momentum, corrections to Eq. ^ arise from the fact that electrons in matter are 
bound to atoms. This implies that there is transverse momentum in the collision (Bloch 
corrections) and that the electron may have momentum comparable to the incident particle 
(shell corrections). These corrections are significant only when the momentum of the incom- 
ing particle is comparable to the electron momentum. As discussed below, we will neglect 
contributions to the slepton range from very low velocities. For (9(100) GeV sleptons trav- 
eling at the velocities we include, the Bloch and shell corrections may be safely neglected. 
The Barkas effect, which introduces a dependence on the sign of the charge of the incoming 
particle is also significant only for very low (3, and may be safely neglected for the velocities 
we include. 

The Bethe-Bloch equation Eq. ^ with 6 as given in Eq. (jZj) is accurate down to /3 ~ 0.05. 
Below this velocity, experimental data are fit by parameterization schemes. We do not have 
this luxury because, of course, there are as yet no experimental data for sleptons. Although 
we could use the parameterizations adopted for standard model particles, it has been found 
that the heavier the particle the worse these models are at describing the low momentum 
behavior js^ . 

Rather than grapple with how sleptons should behave at /5 < 0.05, we adopt the following 
procedure for determining the slepton range in matter. For a certain low velocity, dE/dx 
in Eq. ^ peaks and then rapidly drops to zero. We denote this velocity /Jpeak! its value is 
typically ~ 0.01. At /?peak, the value of dE/dx is very high, and we expect that the distance 
the slepton travels as it slows from Ppeak to thermal equilibrium is negligible relative to the 
distance it traveled in slowing down from its initial velocity to Ppeak- We further take the 
continuous slowing down approximation. With these assumptions, the range i? in g cm~^ 
for a slepton with energy E' is 

RiE') = / dE—, ^ ^ , (8) 

where 6 is given in Eq. ((7j), and M + 6M is the energy at which the slepton has velocity 
Ppeak- The values used for 6M, I, E^, and hujp for the particular case of a 219 GeV slepton 
are given in Table H] 

The Bethe-Bloch equation gives an average value of the stopping power, but there are 
always fluctuations about these values. For thin materials the most probable energy loss 
may be very different from the mean energy loss [sj . However, if the path length in material 
is large enough, the energy loss distribution is Gaussian. The criterion for Gaussianity has 
been found in Ref. [s^ to be 

•^-^^^^^ >10, (9) 



where 



T 

T = , ^ , (10) 



6 



0.2 0.5 1.0 2.0 5.0 10.0 220 230 240 250 

p/M E (GeV) 



FIG. 1: The range to mass ratio R/M as a function of p/M = (3^ (left), and, for the specific case 
of a slepton with mass 219 GeV, the range as a function of energy (right). Results are given for 
water (solid) and lead (dashed). 

and X is the path length. For LHC and ILC detectors, n ^ (9(100). Given that the energy 
distribution of produced sleptons is not wildly fluctuating, inclusion of this Gaussian range 
distribution will have a negligible impact on our results, and we do not include it. 

Using the procedure described above, the range to mass ratios R/M for lead and water 
are given in Fig. [T] These results agree beautifully with published results jHf . The range as 
a function of energy for the particular case of a slepton with mass 219 GeV is also given. A 
slepton with this mass and energy 240 GeV (250 GeV) travels about 10 meters (21 meters) 
in water before stopping. 

III. SLEPTON TRAP OPTIMIZATION 

We will consider the possibility of trapping sleptons in water placed just outside a collider 
detector. Of course, any material may be used, and our analysis, at the level we have pursued 
it, is valid for any material. We consider water to be promising, however, as it is potentially 
feasible to increase the concentration of sleptons and/or move it to a place where slepton 
decays may be observed in a background-free environment. 

We would like to optimize the placement and shape of the water tank. If the gravitino 
LSP scenario is realized in nature, given the implications described in Sec. IVH we consider 
the importance of slepton trapping studies to be sufficient to enlarge detector halls. We 
therefore do not consider constraints from existing detector halls in considering trap geome- 
tries. Even with this simplification, however, a detailed discussion of optimization requires 
careful accounting of the various LHC and ILC detector component geometries, costs, and 
other factors. 

Rather than undertake such a detailed study, we consider here a simple detector model 
to highlight the physics that may be explored in more detail in following studies. We 
characterize the properties of the inner detector by two parameters: r^n, the distance from 
the interaction point (IP) to the outside of the detector, and rf^ = Jpdl, the density- 
weighted distance between the IP and the outside of the detector, typically measured in 
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Interaction Point 



FIG. 2: Diagram of the slepton trap geometry. The trap is assumed to be a spherical shell with 
inner radius r-m, and depth d as shown. The angular parameters ^ A (cos 9) and Acf) of Eq. ((TT|) are 
also indicated. 

meters water equivalent (mwe). We will assume that Tin and are independent of polar 
angle 6 and azimuthal angle 0; that is, we model the inner detector as spherically symmetric 
in both size and material depth. 

The amount of energy lost by a slepton traveling through a detector is determined by 
rj^. (For realistic detector sizes, the tracks of sleptons with sufficient energy to pass through 
the detector have negligible curvature.) The radius of the LHC detectors is approximately 
12A/ in the direction perpendicular to the beam line, where A/ is the nuclear interaction 
length js^ [s^. This number depends on rapidity, but we take this minimal value for 
our spherical detector. The ILC detector is expected to be slightly smaller, but sheets of 
lead or other material can always be used to increase the effective radius to that of the 
LHC detectors. As we will see, this may be advantageous. Given these considerations, we 
assume for this study that the energy loss at the LHC and ILC is approximately that of a 
particle traveling through 12A7- of water. Since Xj = 83.6 g cm~^ for liquid water, this sets 
r'^ = 10 mwe. 

The number of trapped sleptons is, of course, maximized by placing the water tank 
as close to the IP as possible.^ In addition, as we will see, the polar angle distribution 
of slow sleptons, that is, those that have a chance of being trapped, is either uniform or 
peaked perpendicular to the beam line at polar angle 6 = 7i/2. To optimize the water tank 
placement, then, we consider the family of tank geometries specified by 



where (r, 6, 0) are polar coordinates centered on the IP, and d is the tank's radial depth. 



This assumes that no additional material is added between the detector and the trap. It may be advanta- 
geous to place the trap farther away if a material that slows sleptons is added between the detector and 
the trap. An example is discussed in Sec. 



rin < r < Tin + d 




(11) 
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FIG. 3: The depth d in meters of a 1 kton water trap in the (A(cos 0), A(/)) plane for = 10 m. 



This geometry is illustrated in Fig. El All sleptons with range less than r"^^ are trapped 
in the detector. However, all sleptons with range between rj'^ and rf^ + d, and polar and 
azimuthal angles in the ranges given in Eq. (fTT|) are caught in the water tank.^ The water 
tank volume is 

V=\ [(rin + df - rfj A(cos^)A<^ . (12) 

In summary, the number of trapped sleptons is 

iV(r,A(cos0),A0,ri„O , (13) 

where the first four parameters determine the depth d through Eq. ()12j). Sample depths for 
a 1 kton trap {V = 1000 m^we) are shown in Fig. El In the following sections, we fix r^^ and 
to appropriate values and choose three representative sizes V. We then scan over all 
possible values of A(cos^) and A0 to maximize N. In this way, we determine the optimal 
shape for the water tank and the maximal number of sleptons that may be trapped. 



IV. SLEPTON TRAPPING AT THE LHC 

The LHC is scheduled to collide protons with protons at ^/s = 14 TeV beginning in 
2007-08. Its initial luminosity is expected to be 10 fb~^/yr, growing to 100 fb~^/yr. We 
present results below for 100 fb~^, a high luminosity year. 

The ATLAS and CMS detectors are cylindrical with radii 10 m and 8 m, respectively. 
As explained in Sec. 11111 we choose r^^ = 10 m and rj^ = 10 mwe. Although we consider 
only ATLAS and CMS, other detectors could provide promising opportunities. For example, 
LHCb is an asymmetric detector and allows for smaller values of and r^^. This may make 
improved results possible and is a possibility well worth considering. 

In the gravitino LSP scenario, all supersymmetry events produce two long-lived NLSP 
sleptons. The dominant source of NLSPs at hadron colliders is typically pair production 



^ In practice, for the Monte-Carlo simulations below, we smooth distributions in (f> by including all events 
that pass the r and cos 9 cuts with weight A(j>/2Tr. 
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FIG. 4: Representative superpartner masses as a function of in minimal supergravity with 

fixed mo = 0, Aq = 0, tan/3 = 10, and fj, > 0. The supersymmetry parameter n, which governs the 
Higgsino masses, is also shown. 



of strongly interacting superpartners. The number of trapped sleptons is therefore model- 
dependent in a complicated way, as it depends not only on the slepton mass but also sensi- 
tively on the masses of colored superpartners and their cascade decay patterns. 
Here we consider minimal supergravity with the following parameters: 

mo = , Ml/2 = 300 - 900 GeV , Aq = , tan/3 = 10 , fi > . (14) 

When the gravitino is not the LSP, this is in the excluded "stau LSP" region. In the present 
scenario with a gravitino LSP, however, these models are allowed, and this one-dimensional 
family of models provides a simple set with which we can explore the prospects for trapping 
sleptons at the LHC. The lower bound on M1/2 is determined by the requirement of a stau 
NLSP. The number of trapped staus rapidly diminishes as M1/2 increases, and the upper 
bound on M1/2 is roughly where only a few staus may be trapped per year. The superpartner 
spectra for various M1/2 in these models are given in Fig. EJ 

The mass spectra of Fig. |3] are calculated by ISASUSY 7.69 with top mass = 
175 GeV [3^. We have also used this package to generate LHC events. 100,000 non- 
standard model events were generated for minimal supergravity models with the parameters 
given above and M1/2 varying from 300 GeV to 900 GeV in 100 GeV increments. Helic- 
ity correlations are not included in ISAJET. However, there are typically several steps in 
decay chains leading to the NLSP, and so we do not expect helicity correlations to have a 
significant impact on the NLSP distributions or on our final results. 

In Fig. we show the energy distribution of NLSP staus for the M1/2 = 600 GeV. 
Although many staus are produced, most of these are extremely energetic and impossible to 
stop in a reasonable distance. Of course, a trap could be set up far from the IP so that the 
intervening earth slows down the slepton, but such a trap would be too far away to have a 
reasonable solid angle coverage for any realistic volume V. 

The cos 6* distribution at M1/2 = 600 GeV is shown in Fig. IHl Given all produced NLSPs, 
the distribution is strongly peaked in the beam directions. However, imposing a cut on 
energy, we find that the slow staus are produced roughly isotropically. Thus, for those staus 
that we might reasonably hope to be trapped, any polar angle is as good as any other. Given 
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E (GeV) 

FIG. 5: The energy distribution of NLSP staus produced at the LHC for integrated luminosity 
100 fb^^ and minimal supergravity with niQ = 0, M1/2 = 600 GeV, ^0 = 0, tan/? = 10, and /u > 0. 
The NLSP stau mass is 219 GeV. 




COS0 



FIG. 6: The cos^ distribution of NLSP staus produced at the LHC for integrated luminosity 
100 fb^^ and minimal supergravity with mo = 0, M1/2 = 600 GeV, ^0 = 0, tan/3 = 10, > 0. The 
NLSP stau mass is 219 GeV. Distributions are given for all staus, Ef < 290 GeV {R < 100 mwe), 
and Ef < 240 GeV {R < 10 mwe). The total distribution is strongly peaked along the beam 
directions, but slow sleptons are produced isotropically. 

the reality that LHC detectors are cylindrical, however, the closest a trap may be placed to 
the IP is at cos 6* = 0, justifying our choice of centering our trap geometries at cos^ = 0, as 
parameterized in Eq. (jllll . 

The number of trapped sleptons for optimized trap shape and placement and various trap 
volumes is given in Fig. [7| The trap is optimized as described in Sec. Illll we scan over all 
possible A(cos^) and A(f), and find the combination that maximizes the number of sleptons 
that stop in the trap. 

We find that only a small fraction of produced sleptons can be trapped. For example, 
for Ml/2 = 600 GeV, 4.2 x 10^ NLSP sleptons are produced, but only 260, 40, and 9 are 
trapped in water tanks of size 10, 1, and 0.1 kton, respectively. For all models considered. 
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100 150 200 250 300 
I^NLSP 

FIG. 7: The number of sleptons trapped per year at the LHC in water tanks of size 10 kton (sohd), 
1 kton (dot-dashed), and 0.1 kton (dashed). The total number of sleptons produced is also shown 
(upper dotted), along with the number of sleptons trapped in the LHC detector (lower dotted). 
The water tank shape and placement have been optimized as described in the text. These results 
assume luminosity 100 fb~^/ yr, Tin = 10 m, = 10 mwe, and minimal supergravity models with 
Ml/2 = 300, 400, . . . , 900 GeV, mo = 0, = 0, tan (3 = 10, and ;U > 0. 

the 10 kton trap is optimized for A(cos6') = A^/vr = 2; for this large volume, the trap 
is sufficiently deep that there is little be gained by making the trap deeper at the expense 
of solid angle coverage. For V = 1 kton, the number of trapped sleptons is optimized by 
traps with less-than- maximal angular coverage and depths of 4 m. Note that, because 
slow sleptons are produced isotropically, the number of trapped sleptons depends, to a good 
approximation, on A(cos^) and A0 only through their product A(cos^)A0. 

Despite the low efficiency for trapping, so many NLSPs are produced at the LHC that 
significant numbers of NLSPs may still be trapped. As anticipated, the results are heavily 
dependent on the overall scale of superpartner masses. For V = 10 kton and mNLSP = 
100 — 300 GeV, the number of trapped sleptons varies from O{10^) to (9(10). For the lighter 
sleptons considered, these results imply that sufficient numbers of sleptons may be trapped 
to do precision studies of slepton decay properties. Assuming that slepton decays may be 
observed in a background-free environment, we expect percent level measurements of slepton 
decay widths. 

Note that for all M1/2, a larger number of sleptons range out in our spherically symmetric 
LHC detector than can be trapped in even the 10 kton water trap. For M1/2 = 600 GeV, 
700 are trapped in the LHC detector itself. The LHC detector has a large volume and 
beneffis from the fact that it begins at the IP and so has a large angular coverage without 
sacrificing depth. Unfortunately, it is not clear to what extent these sleptons may be used — 
their decays are out of time, occur away from the IP, and take place in an environment with 
significant cosmic ray background. Given the large number of sleptons that are automatically 
trapped in the LHC detector itself, however, it is certainly worthwhile to explore ways to 
exploit them. 
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V. SLEPTON TRAPPING AT THE ILC 



In its first stage, the ILC will collide electrons and positrons with center of mass energies 
up to 500 GeV. In this first stage, the luminosity has been estimated to be 340 fb~^/yr for 
the TESLA design and 220 fb^Vyr for the NLC/JLC For this study, we assume 
luminosity 300 fb~^/yr. As in the LHC analysis above, we present results for one year of 
running. 

At present ILC detectors are expected to be slightly smaller than their LHC counterparts. 
To be conservative, we assume = 10 m and = 10 mwe, the same parameters we 
assumed in the LHC case. Of course, if the detector is smaller than this, it can always be 
supplemented by adding plates of lead, for example, to mock up these parameters. As we 
will see, in the ILC case, such an approach may in fact enhance our results very significantly. 

At the ILC, scanning over supersymmetry models with a broad range of superpartner 
mass scales, as done in the LHC analysis above, is not particularly informative. Models 
with heavy superpartners are simply out of reach, and no sleptons may be produced, much 
less trapped. On the other hand, for models with superpartners within reach, the ILC beam 
energy may be tuned to optimize the number of trapped sleptons, to some extent offsetting 
variations in the scale of superpartner masses in these models. As we will see, the crucial 
feature is not the exact mass of the slepton NLSP, but rather the presence of other nearly 
degenerate superpartner states. 

For the ILC, then, we limit our analysis to two models. In the first, which we denote 
"NLSP only," the only superpartner within reach of the ILC is an NLSP f/j with mass 219 
GeV. This is representative of the minimal case where the gravitino LSP scenario may be 
probed at the ILC. Of course, in many realistic models, there are a number of other super- 
partners, notably other sleptons, fairly degenerate with the NLSP. We therefore consider 
also a second model, which we denote "mSUGRA," which is minimal supergravity with 
Ml/2 = 600 GeV, Aq = 0, tan/3 = 10, and /i > 0. This model contains not only the 219 GeV 
fpt of the "NLSP only" model, but also right-handed selectrons, right-handed smuons, and 
a neutralino within the kinematic reach of a 500 GeV ILC. The mSUGRA model is one of 
the family of models considered previously in the LHC analysis, allowing us to compare the 
LHC and ILC at one particular model point. Because it contains the "NLSP only" model 
as a subset, it also allows us to see the effect of having other accessible and fairly degenerate 
superpartners. The accessible standard model superpartners of the two models and their 
masses are: 



242.9 GeV 
mg«,m^^ 227.2 GeV 



mSUGRA (15) 



TUfj^ 219.3 GeV } NLSP only 

We generate 10'^ non-standard model ILC events for the mSUGRA model with ISASUSY 
7.69 with nit = 175 GeV j3^. Events for the "NLSP only" model are compiled by selecting 
the prompt stau events from this event sample. We choose beam width 0.12 mm, and 
beamstrahlung parameter T = 0.1072, and allow the subprocess energy to vary over the 
entire range from 2mNLSP to a/s. 

In the "NLSP only" model, the NLSP staus are produced through e~^e~ — 7, Z — 
f+f". The stau energy distribution is therefore given by the beam energy modified by 
initial state radiation (ISR) and beamstrahlung. An example with ^/s = 500 GeV is given 
in Fig. ISl The stau polar angle distribution is also given in Fig. |H| Despite ISR and 
beamstrahlung, it retains the sin^ 9 shape of the parton-level process. The best place to 
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FIG. 8: Energy (left) and cosO (right) distributions for NLSP staus produced at the ILC assuming 
the "NLSP only" model, in which the only accessible superpartner is the NLSP stau with mass 
219 GeV. Results are for ^/s = 500 GeV and integrated luminosity 300 fb~^. 
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FIG. 9: The number of sleptons trapped per year at the ILC in 10 kton (solid), 1 kton (dot- 
dashed), and 0.1 kton (dashed) water traps. The total number of sleptons produced is also shown 
(upper dotted) along with the number of sleptons trapped in the ILC detector (lower dotted). The 
trap shape and placement have been optimized, and we assume Tin = 10 m and rj^ = 10 mwe, 
luminosity 300 fb~^/yr and the NLSP only model, where the only accessible superpartner is a 219 
GeV NLSP stau. 



trap sleptons is therefore perpendicular to the beam line, justifying our choice of centering 
our trap geometries at cos^ = 0, as parameterized in Eq. (fTT|) . 

The number of trapped sleptons for various trap sizes as a function of center-of-mass 
energy ^/s is given in Fig. IHl For -y/i < 475 GeV, no staus escape the ILC detector. At 
^/s = 475 GeV, however, sleptons in the sharp peak of the energy distribution escape the 
ILC detector and may be caught in a fairly thin water tank placed just outside the ILC 
detector. 

As evident in Fig. El for the 1 and 0.1 kton water traps, the number of trapped staus is 
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FIG. 10: The number of trapped staus at the ILC in the (A (cos 6), A(j)) plane for the "NLSP only" 
model, ^/s = 475 GeV, and integrated luminosity 300 fb~^/yr. 



maximized when the beam energy is tuned to produce staus that just barely emerge from the 
ILC detector. The dependence on trap parameters is illustrated in Fig. [TUl The optimized 
trap configuration has A(cos^) ~ 1 and A0 = 2tt; because the stau distribution is peaked 
at cos 6* = 0, it is beneficial to sacrifice coverage at high rapidity to make the trap deeper. 
For a 10 kton trap, the trap is sufficiently thick that the best results are achieved for slightly 
higher beam energies where more of the ISR/beamstrahlung tail may be caught. For a 10 
kton trap, we find that the optimal trap configuration has A(cos0) ~ 2 and A0 = 27i. 

For all trap sizes, however, the number of trapped sleptons is maximized for beam energies 
near the 475 GeV threshold. The tunable beam energy and well-defined initial state are well- 
known virtues of the ILC, but these are exploited in a qualitatively novel way here to produce 
slow NLSP sleptons that may be easily caught. Relative to the case at the LHC, a much 
larger fraction of the produced staus can be caught. For example, of the 2650 staus produced 
at ^/s = 482 GeV in a year, 2000 staus may be trapped in an optimized 10 kton water trap! 
Such results imply promising prospects for slepton trapping even in the minimal case when 
the only superpartner accessible at the ILC is the NLSP slepton. 

Before considering the mSUGRA model, we note that the number of produced staus 
continues to rise well beyond its value at = 475 GeV. This suggests that our results may 
be improved significantly by placing some dense material between the ILC detector and the 
water tank. By adding material depth to the ILC detector, the threshold at which sleptons 
just barely emerge is moved to higher y/s where the stau pair production cross section is 
higher. For a dense material, such as lead, this can be achieved without increasing r^^ much. 
Such a strategy may in any case be required to smooth out variations in inherent in 
realistic detectors. Although we have not investigated this in detail, we expect that a large 
enhancement may be possible. 

We now turn to the mSUGRA model. In this model, NLSP staus may again be 
produced directly, but now they may also be produced in several other ways: first, by 
e~^e~ —>■ e~ , fi'^ fi~ followed by e — *• err and ft — > /irr, and second, by e+e" — > XX: fol- 
lowed hy X ^ T'r, or through the cascade x ^ e, /i — >^ f. The energy and cos^ distributions 
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FIG. 11: Energy (left) and cos^ (right) distributions for NLSP staus produced at the ILC in the 
mSUGRA model with mo = 0, M1/2 = 600 GeV, Aq = 0, tan/3 = 10, > 0. Results are for 
y/s = 500 GeV and integrated luminosity 300 fb^^/yr. 

of NLSP staus in the mSUGRA model are shown in Fig. As evident in Fig. ^2 the 
additional sources of staus increase the total number of staus significantly, but just as signif- 
icant, the cascade decays produce a broad and flat tail in the energy distribution extending 
nearly down to mf. The cos^ distribution is nevertheless still peaked at cos 6* = 0. 

The number of trapped staus per year for the mSUGRA model are given in Fig. El The 
presence of additional accessible superpartner states has a significant impact — for all trap 
sizes considered, large numbers of staus may be trapped even for beam energies well above 
475 GeV. This is a consequence of the broad energy distribution of NLSP staus, which in 
turn follows from the existence of other fairly degenerate superpartners. 

In the mSUGRA case, we may also compare these ILC results directly with results from 
the LHC analysis given above. The LHC results for M1/2 = 600 GeV from Sec. llVl are given 
by the lines marked "LHC" in Fig. ^1 Comparing results of similarly-sized water traps, we 
find that for this particular model, the ILC will be able to trap a factor of ~ 10 more staus 
than the LHC. Again, for the reasons discussed above, a significant enhancement of these 
ILC results may be possible if one considers inserting lead between the ILC detector and the 
water trap, and as many as C(10'') staus may be trapped by reasonably sized water traps 
using this more general approach. 

VI. IMPLICATIONS AND CONCLUSIONS 

Although the implications of supergravity for cosmology and particle physics have been 
considered in great detail for decades, most work has been centered on scenarios in which 
the LSP is a standard model superpartner. Here we have explored the gravitino LSP sce- 
nario. Recent work has found significant cosmological motivations for this possibility, as 



Helicity correlations between production and decay are not included in our event generation. These are, 
of course, absent for scalar particles, but may modify both the energy and cos 9 distributions for staus 
produced through x P&ir production. 
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FIG. 12: The number of sleptons trapped per year at the ILC in 10 kton (solid), 1 kton (dot-dashed), 
and 0.1 kton (dashed) water traps. The total number of sleptons produced is also shown (upper 
dotted) along with the number of sleptons trapped in the ILC detector (lower dotted). The trap 
shape and placement have been optimized, and W6 cLSSume t\yy — 10 m and. — 10 mwG, luminosity 
300 fb-i/yr and the mSUGRA model with mo = 0, M1/2 = 600 GeV, ^0 = 0, tan/3 = 10, ^ > 0. 
For comparison, the number of trapped staus at the LHC for various trap volumes is also shown. 



the gravitino may explain dark matter, and the scenario may resolve current difficulties in 
Big Bang nucleosynthesis and with leptogenesis. We considered here a novel implication for 
collider physics, namely that NLSP sleptons may be collected in water traps before their 
decays to the gravitino. These sleptons may then be concentrated and transported to some 
quiet environment where their decays may be studied in detail. 

By optimizing the water trap shape and placement and considering a variety of sizes, we 
have ffist explored the prospects for trapping sleptons at the LHC. The number that may 
be trapped is highly model-dependent. For minimal supergravity with mg = 0, we find that 
as many as 10^ staus may be stopped in a 10 kton trap when the sleptons have mass around 
100 GeV. This is as light as is allowed by current bounds. For a less optimistic scenario, say, 
with 219 GeV staus, hundreds and tens of sleptons may be trapped each year in 10 kton 
and 1 kton traps, respectively. 

These results may be improved significantly if long-lived NLSP sleptons are kinematically 
accessible at the ILC. For the identical case with 219 GeV sleptons discussed above, 0(1000) 
sleptons may be trapped each year in a 10 kton trap. If only the NLSP is accessible, this 
result may be achieved by tuning the beam energy so that produced NLSPs barely escape 
the ILC detector. The ability to prepare initial states with well-known energies and the 
flexibility to tune this energy are well-known advantages of the ILC. Here, these features are 
exploited in a qualitatively new way to produce slow sleptons that are easily captured. 

If there are additional superpartner states accessible at the ILC, even tuning the beam 
energy is not necessary. The cascade decays of other superpartner states produce a broad 
distribution of slepton energies, and so for a broad range of beam energies, some sleptons 
will be captured in the trap. We have noted also that, by considering the slightly more 
general possibility of placing lead or other dense material between the ILC detector and the 
slepton trap, an order of magnitude enhancement may be possible, allowing up to (9(10'^) 
sleptons to be trapped per ILC year. 
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The analysis here is vahd with minor revisions for traps composed of any material. For 
concreteness, however, we have considered traps composed of water tanks, with the ex- 
pectation that sleptons caught in water will be easily concentrated and/or moved to quiet 
environments. 

These results imply that high precision studies of slepton decays may be possible. There 
are many significant implications of such studies. These have been considered in detail 
in Refs. j^, [l^. Briefly, simply by counting the number of slepton decays as a function 
of time, the slepton lifetime may be determined with high accuracy. Given thousands of 
sleptons, we expect a determination at the few percent level. The slepton decay width of 
Eq. (jSj) is a simple function of the slepton and gravitino masses, and the slepton mass will 
be constrained by analysis of the collider event kinematics, a measurement of the slepton 
width therefore implies a high precision measurement of the gravitino mass and, through 
Eq. (fP), the supersymmetry breaking scale F. Such measurements will provide precision 
determinations of the relic density of superWIMP gravitino dark matter, the contribution 
of supersymmetry breaking to vacuum energy, and the opportunity for laboratory studies 
of late decay phenomena relevant for Big Bang nucleosynthesis and the cosmic microwave 
background. 

The gravitino mass may also be determined, although not necessarily on an event-by- 
event basis, by measuring the energy of slepton decay products. This provides a consistency 
check. Alternatively, these two methods, when combined, determine not only itlq, but also 
the Planck mass M*. This then provides a precision measurement of Newton's constant on 
unprecedentedly small scales, and the opportunity for a quantitative test of supergravity 
relations. 
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Note added: E-print hep-ph/0409248 [39|, in which long-lived sleptons at colliders are also 
studied, appeared as this work was being finalized. This paper studies the possibility of 
using an active detector and also considers the idea of using e~e~ collisions at the ILC. 
The possibility of transporting trapped sleptons to a low background environment, and the 
material-independent analysis of optimizing trap shape and placement, both discussed at 
length here, were not addressed. 
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